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We study noncentrosymmetric effects on resonant x-ray scattering (RXS) in magnetite. 
The noncentrosymmetry at A sites in spinel structure makes the Ap states strongly hybridize 
with the 3d states through neighboring oxygen 2p states, giving rise to the non-vanishing 
contribution of the dipole-quadrupole (E1-E2) process in the RXS spectra. We substantiate 
this observation by introducing a microscopic model of a FeC>4 cluster with multiplets and 
the Ap band. We show that the hybridization changes its sign between two kinds of A sites 
and accordingly the local amplitude from the E1-E2 process changes its sign, resulting in 
non-vanishing RXS intensities at forbidden spots (002) and (006) in the pre-edge region in 
agreement with the experiment. A large dependence of the pre-edge intensity on the direction 
of the applied magnetic field is predicted as a consequence of breaking both centrosymmetry 
and time-reversal symmetry. Furthermore we analyze the intensity difference between two 
opposite directions of the applied magnetic field at the allowed spot (222) in connection with 
the experiment. We obtain the intensity difference of a "dispersion" form, which resembles the 
observed spectra at the Mn pre-if-edge in MnCr204 but is quite different from the observed 
one in magnetite. Although the observed spectra are claimed to arise from "magnetoelectric" 
amplitude, we argue that this claim has no ground. 

KEYWORDS: resonant x-ray scattering, magnetite, local noncentrosymmetricity, E1-E2 tran- 
sition, magnetoelectric effect, XMCD, MnCr2C>4 



1. Introduction 

Resonant x-ray scattering (RXS) has been widely used to investigate different kinds of 
orders, such as charge, magnetic and orbital orders in crystals, 1-5 since the strong resonance 
makes the signal sensitive to the ordered structure. The K-edge resonance is usually used in 
transition metals in order to observe signals at superlattice spots associated with the order 
parameter. This is because the corresponding x-ray wavelength matches the period of long 
range orders, which is usually an order of atomic distance except for the case of long period. 

The RXS amplitude is given by a sum of atomic amplitudes with appropriate phases. 
Each atomic amplitude is described by a second order process. One of the most dominant 
processes in transition-metal compounds is the dipole-dipole (El-El) process that the Is 
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electron is excited to the Ap states by absorbing x-ray and then the Ap electron is recombined 
with the ls-core hole by emitting x-ray. Since the Ap states are extended in space, they are 
easily influenced by the electronic structure at neighbors to the core-hole site; the Ap states 
are modulated by the lattice distortion through the hybridization to neighboring oxygens, 
giving rise to the signal at superlattice spots. Therefore, the RXS signal at the superlattice 
spots arises from the variation of the Ap states in accordance with the long-range order, and 
provides an indirect proof of the order which is usually constructed by 3d states. Such a view 
has been confirmed by theoretical analyses 6-8 in connection with the RXS experiment for 
LaMn0 3 . 3 

There appear sometimes extra signals with energy below the K-edge, called the pre-edge 
signals, which will be mainly discussed in the present paper. Since the pre-edge energy is 
close to the energy exciting an electron from the Is state to the 3d states, the signal could 
be naturally interpreted as arising from the quadrupole-quadrupole {E2-E2) process in which 
the Is electron is excited to the 3d states by absorbing x-ray and then one of 3d electrons 
is combined with the core hole by emitting x-ray. However, the pre-edge signal could also 
be generated from the El-El process, since the p-symmetric states with respect to the core- 
hole site can be constructed from 3d states at neighboring transition-metal atoms. 6 ' 9 These 
two origins may be distinguished by different peak positions, that is, the peak in the E2-E2 
process is expected to be located at the region around several eVs lower than that in the 
El-El process, since the relevant 3d states in the E2-E2 process is on the core-hole site, and 
is strongly attracted by the core-hole potential. 

The situation may become quite different when the centrosymmetry is locally broken. In 
such circumstances, the Ap states could hybridize with the 3d states on the same site through 
the hybridization to neighboring oxygen 2p states, and thereby the dipole-quadrupole (E1-E2) 
process could contribute to the pre-edge signals. Such a presence of the E1-E2 process has 
been recognized by the experiment of K2Cr04 10 and by the numerical calculation for Ge. 11 
It has been much debated for V2O3. 12-15 Furthermore, for magnetic materials, the atomic 
amplitude of the E1-E2 process could depend on the direction of the local magnetic moment 
due to the spin-orbit interaction (SOI), and thereby the pre-edge signals could depend on the 
direction of magnetic moment. Since the direction of magnetic moment could be controlled 
by applying the external magnetic field, we could observe such a dependence by changing 
the external magnetic field. Actually, such signals have been observed 16 and analyzed in a 
multiferroic system GaFeOs, 17,18 and also have been discussed in other situations. 19,20 

The pre-edge signals have also been observed in magnetite, Fes04, at forbidden spots 
of scattering vectors (002) and (006). 21,22 In addition, the intensity difference with changing 
direction of the external magnetic field has been measured at spots (222), (333) and (444). 23,24 
The purpose of this paper is to analyze the pre-edge signals in magnetite through a quantitative 
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calculation of the spectra and to elucidate the mechanism from a microscopic viewpoint. 

The magnetite is the first magnetic material known to the mankind. Its crystal structure is 
the inverse spinel, consisting of iron sites tetrahedrally surrounded by four oxygens (A sites) 
and those surrounded octahedrally by six oxygens (B sites), as shown in Fig. 1. Since the 
centrosymmetry is locally broken at the tetrahedral sites, those pre-edge signals are thought 
to be related to breaking both centrosymmetry and time-reversal symmetry. Analyses based 
on the microscopic electronic structure, however, have not been worked out yet. We construct a 
definite model that the 4p states form an energy band with wide width and hybridize strongly 
with the 3d states through neighboring oxygen 2p states. The multiplet structures are taken 
into account in the 3d 5 - and 3<i 6 -configurations. Applying the resolvent formalism 25 to this 
model, we calculate the local electronic structure around the tetrahedron sites and thereby 
the atomic amplitudes of RXS. 




It is important to recognize that there are two kinds of tetrahedron sites denoted as Ai 
and A2, that is, one is transformed into the other by space inversion with respect to the center 
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of the tetrahedron (see Fig. 2). We find that the effective hybridization between the 4p and 
3d states via oxygen 2p states changes its sign between the Ai and A2 sites, leading to a sign 
change in the atomic amplitude of the E1-E2 process. This is a key point to explain how the 
pre-edge signals come out. At spots (002) and (006), the contributions from the El-El and 
E2-E2 processes as well as Thomson scattering are canceled out in the a — a 1 channel, and that 
of the E1-E2 process at the A sites only survives in the total scattering amplitude. 26 We obtain 
the pre-edge spectra as a single peak as a function of the photon energy, in agreement with the 
experiment. 21 ' 22 Furthermore, we calculate explicitly the dependence on the direction of local 
magnetic moment in the atomic amplitude. The depending parts are found about an order of 
magnitude smaller than non-dependent ones. From this calculation, we obtain the intensity 
differences at spots (002) and (006) when the direction of the magnetization is changed from 
the [1,-1,0] direction to the reverse, which shape looks like an "absorption" type as a function 
of photon energy. It may not be hard to detect these signals, since the magnitudes are about 
1/5 to the corresponding pre-edge intensity peaks. 




Fig. 2. Two types of tetrahedrons, A\ (left) and A 2 (right). 



We also analyze the dependence on the external magnetic field at spots (222), (333) and 
(444) in connection with the recent experiment. 24 These are allowed spots, where the Thomson 
scattering amplitude is dominant. Focusing on the E1-E2 process at the A sites, we calcu- 
late the intensity difference between two opposite directions of the applied magnetic field. 
The difference arises from the interference between the Thomson scattering amplitude and 
the E1-E2 amplitude. We show that the intensity differences are nearly the same magnitude 
at both (222) and (333) spots but no difference at (444), and that the shapes as a function 
of photon energy look like a "dispersion" form concentrated in the pre-edge region. In the 
experiment by Matsubara et al., 24 however, the intensity difference at the (222) spot is dis- 
tributed over the region much wider than the pre-edge region and two orders of magnitude 
larger than that at (333) and the calculated values. Also, the spectral shape is quite different 
from the "dispersion" form. Matsubara et al. claimed that the difference at (222) arises from 
a "magnetoelectric" amplitude, that is, a consequence of breaking both centrosymmetry and 
time-reversal symmetry. We argue that this claim has no ground. Our finding of a "dispersion" 
form for the intensity difference has been observed in the experiment at the Mn pre-i^-edge 
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in MnCr2 04, 27 where Mn atoms are occupying at the A sites in spinel structure. Since the 
pre--KT-edge signal selects the contribution from the A sites, this experiment suggests that the 
calculated spectra correspond to the signal from the E1-E2 process. 

This paper is organized as follows. In Sec. II, we briefly summarize fundamentals of mag- 
netite. In Sec. Ill, we introduce the model Hamiltonian and study the electronic structure 
around the A sites. In Sec. IV, we describe the excited states involving a ls-core hole by 
employing the resolvent formalism. In Sec. V, we calculate the absorption spectra and discuss 
the x-ray magnetic circular dichroic (XMCD) spectra. We calculate the RXS spectra in com- 
parison with experiments. The last section is devoted to concluding remarks. The geometrical 
factors are summarized in Appendix. 

2. Fundamentals of magnetite 

The crystal structure of magnetite is the inverse spinel with the lattice constant a® = 
8.396 A, as shown in Fig. 1. The unit cell contains 24 iron atoms; 8 iron atoms are tetrahedrally 
coordinated by 4 oxygens (A sites) and 16 iron atoms are octahedrally coordinated by 6 
oxygens (B sites). Note that two types of tetrahedrons exist within the A sites (Ai and A2), 
as shown in Fig. 2. It will be clarified in subsequent sections that the subtle difference in 
electronic structures between the Ai and A2 sites leads to important consequences on the 
RXS spectra. With disregarding small crystal distortion, we have the position vectors of Fe 
atoms in the unit cell, 
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(2.1) 

Iron atoms at the A sites are nominally Fe 3+ , while those at the B sites are a 1 : 1 
mixture of Fe 3+ and Fe 2+ . Note that Fe 3+ atoms are in the 3d 5 -configuration with the spin 
angular momentum S = 5/2, and that Fe 2+ atoms are in the 3<i 6 -configuration with 5 = 2, 
according to the Hund rule. The hybridization between the 3d states and oxygen 2p states may 
change but slightly the situation. The local magnetic moments are ferromagnetically aligned 
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within the individual A and B sites, while those at the A and B sites are antiferromagnetically 
aligned. As a result, a net magnetization remains finite, that is, the magnetite is a ferrimagnet 
at low temperatures. The Curie temperature is as high as around 850 K. In addition, a metal- 
insulator transition so called Verway 28 transition takes place around T = 120 K. This may be 
related to charge and orbital orders, which study is outside our scope. 29-31 

3. Electronic Structures around the A sites 

In this section we focus on the electronic structure around the A sites which have no 
centrosymmetry. In particular, we are interested in the excited states having one Is core hole 
and one 4p electron in accordance with the El process and those having one Is core hole and 
one more electron in the 3d states in accordance with the E2 process. 

3. 1 Crystal electric field 

We start by examining the crystal electric field (CEF) to look for noncentrosymmetric 

effects. Let charge q be placed at the apexes of a tetrahedron. The electrostatic potential 

4>(x, y, z) is expanded around the center as 

4q 20 q 35qf 4 4 4 3 4 

4> (x,y,z = — T -7=— xyz - — — [x +y +z - -r 

where r = \J x 2 + y 2 + z 2 with ro being the distance between the origin and the apexes. 
The last term is well known to represent a split of the energy level of 3d states. The second 
term, which is usually neglected, gives rise to a hybridization between 3d and 4p states. This 
coupling comes out because of noncentrosymmetry, but it is much smaller than the same type 
of coupling arising from the hybridization between the 3d and oxygen 2p states and between 
the 4p and 2p states. The sign — (+) of the coupling is taken for the A\(A2) sites. 

3.2 Effective hybridization between the 4p and 3d states 

Now we discuss how the 4p states could hybridize with the 3d states in the absence of 
centrosymmetry. Let 3d wavefunctions be ip^ 2 , i>1%_ r 2-, Y^f, V'zrc) an d V'xj/j an< ^ wavefunc- 
tions be ipt p , ipy P , and ipi p . They are all real and normalized, and have symmetries described in 
the subscript. Each state could hybridize with a state constructed from a linear combination 
of oxygen 2p wavefunctions at apexes. These oxygen wavefunctions have the same symmetry 
as their partner of hybridization, which are denoted as ip^_ y2 , ^ z 2_ r 2i tpyz, ip'zx, ipxy, '4'x P , 
ipy P , and ipz P - Using the Slater-Koster two-center integrals given in Table I, we evaluate the 
strength of hybridization between the 3d and 2p states and between the 4p and 2p states, 

tT 2p = (^l<b 2p l«_ y2 > 

= «-, 2 |<b 2P !V'J 2 - r2 ) = l-34eV (3.1) 
= (^l<b 2p |^)=T2.33eV, for A ± (A 2 ), (3.2) 
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= (^ P \H^ 2 V/) = -4.36 eV, (3.3) 

where H^ 2p and H^ 2p are the hybridization energies between the 3d and 2p states and 
between 4p and 2p states, respectively. The sign — (+) in eq. (3.2) corresponds to the Ai (A2) 
sites. Note that ipy% and ip 2p are not identical but have a finite overlap. The same is true for 
ip 2p and ipl P and for ip 2p and ipz P , respectively. The overlap is evaluated as 

S = <V#Wg> = = (€ P I^) = "0.748. (3.4) 

Needless to say, these values depend on the phase of wavefunctions constructed from oxygen 
2p orbitals, but the effective hybridization between the Ap and 3d states are independent of 
the phase, since it is proportional to (ipi P \H^ 2p \^x P )(^x P \^yz)(^yz\H^ 2p \^yt)- Its si S n is 
opposite between the A\ and A2 sites. This corresponds to the sign change of the second term 
of eq. (3.1) in the CEF. 
3.3 Hamiltonian for a FeO^ cluster 

Now that the 4p states could hybridize with the 3d states through oxygen 2p states, 
we include oxygen states into our model, in addition to the Is, 3d, and Ap states, in order 
to describe the electronic structure around the A sites. For this reason, we consider the 
Hamiltonian of a Fe04 cluster at the A sites, 

H = H 3d + H 2p + H™~ 2p + H ls + H 4p + H* p ~ 2p , (3.5) 

where 

H 3d = ^Ef n S ma d ma + - ^2 9{v\V2\v-iv±)dl^d Vl ,d V2> 

ma v\V2 V4 

+ (3d ( ma \ l ' s \ m ' a ') d ma d m>a> + (H xe + H ext ) ■ ^ (s^/G^G^/ , (3.6) 

mm' era' maa' 

H 2p = J2E p pl aPma , (3.7) 

ma 

<b 2P = E^" 2P ^m. + H.c, (3.8) 

ma 

H ls = e ls J24s*, (3.9) 

a 

H* P = Y.^)p'laAr,a, (3-10) 
kr/<7 

rja 

The H sd describes the energy of 3d electrons, where d ma represents an annihilation opera- 
tor of a 3d electron with spin a and symmetry m (= x 2 — y 2 , 3z 2 — r 2 ,yz, zx, xy) at the center. 
The 3d energy level is split by the CEF energy WDq. The second term in eq. (3.6) repre- 
sents the intra-atomic Coulomb interaction with the interaction matrix element g {y\v<i\ 
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in terms of F°, F 2 , and F 4 (v stands for (m, a)). The third term in eq. (3.6) represents the 
SOI for 3d electrons with the SOI coupling £3^. We evaluate atomic values of F 2 , F 4 , and C3d 
using the wavefunctions within the Hartree-Fock (HF) approximation, 32 and multiply 0.8 to 
these atomic values with taking account of the slight screening effect. 33 On the other hand, we 
multiply 0.25 to the atomic value for F°, since it is known that F° is considerably screened by 
solid-state effects. The last term in eq. (3.6) describes the energy due to the exchange interac- 
tion from neighboring Fe atoms and the Zeeman energy with the external field, where (s) aa i 
represents the matrix element of the spin operator of the 3d electrons. The exchange field H xc 
here has a dimension of energy, and is an order k^T c with T c = 850 K. The external field 
H cxt is assumed to be much smaller than H xc but to be larger than the magnetic anisotropy 
energy. Therefore it has a role to align the magnetization to the field. 

The H 2p represents the energy of oxygen 2p electrons, where p ma is the annihilation 
operator of the state \tpm) with spin a. The Coulomb interaction is neglected in oxygen 2p 
states. The H^^ 2p represent the mixing energy between the 3d and 2p states, where tm~ 2p 
is the matrix element given by eqs. (3.1) and (3.2). The energy of the 2p level relative to 
the 3d levels is determined from the charge-transfer energy A defined by A = E d - E p + 
15U(3d G ) — 10U(3d 5 ) with E d being an average of E^. Here U(3d 6 ) and U(3d 5 ) are the 
multiplet-averaged d-d Coulomb interaction in the 3<i 6 and 3d 5 configurations, which are 
defined by U = F° - (2/63) F 2 - (2/63) F A with F°, F 2 , and F A . 

The last three terms are added to the Hamiltonian in accordance with the excitation of the 
core electron. The H ls represents the energy of the Is electrons, where s a is an annihilation 
operator of the Is state. The H 4p represents the energy of the Ap states, which form a band 
with energy 64 p (k). The F[f p ^ 2p represents the hybridization between the Ap and oxygen 2p 
states, where p' va is the annihilation operator of Ap electron with symmetry rj = x,y, and z, 
and p' va = (1/y/No) YlkPicna i^o ls t ne number of k). This expression could be changed into 
a form that Ap states hybridize with oxygen states symmetrized as yz, zx and xy: 



with m = yz corresponding to rj = x and so on. Here the matrix element t 4p ~ 2p is renormalized 
as t Ap ~ 2p = t Ap ~ 2p S. We do not explicitly consider the Coulomb interactions between the core 
hole and the Ap and 3d electrons, but we could take account of the main effects by adjusting 
the energy separation between 3d level and e± s , since the Slater integrals responsible to the 
exchange interaction is rather small, G 2 (ls,3<i) = 0.058 eV. 

Table I lists the parameter values used in this paper, which are consistent with the values 
in the previous calculations. 34 




(3.12) 
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Table I. Parameter values for a FeC>4 cluster in the 3d 5 configuration, measured in units of eV. 



F°(3d, 3d) 


6.39 


{pda) 2p ,3d 


-1.9 


F 2 (3d, 3d) 


9.64 


(pdvr) 2 p,3d 


0.82 


F 4 (3d, 3d) 


6.03 


(ppv)2p,4 P 


3.5 


(3d 


0.059 


{PPX)2 P AP 


-1.0 


A 


3.5 


lODq 


-0.7 



3.4 Lowest energy state at the A sites 

Iron atoms at the A sites are nominally the Fe 3+ (3d 5 ) configuration. The hybridization 
between the 3d states and oxygen 2p states makes it mix with the 3d e L configuration, where 
L indicates the presence of a hole in the ligand oxygen orbitals. Preparing 2352 bases in the 
3d 5 + 3d Q L configuration, we represent the Hamiltonian + Hi v + H^^ 2p for H xc along 
the z axis. Diagonalizing numerically the Hamiltonian matrix, we obtain the spin moment 
S = 2.40 and the orbital moment L = 0.0036 in the lowest energy state. These values deviate 
slightly from S = 5/2 and L = in the lowest state of the 3d 5 configuration. Note that these 
values are insensitive to the magnitude and direction of H xc - The weig ht of the 3d 5 and 3d 6 L 
configurations are obtained as 0.795 and 0.205. 

4. Excited states relevant to the K edge RXS 

4 ■ 1 Resolvent formalism 

We use the resolvent formalism in order to describe the excited states containing a Is core 
hole and a 4p electron. It is defined by 

G(z) = [z - H - V]- 1 , (4.1) 

where 

H = H 3d + H 2p + H 3d - 2p + H ls + H 4p , (4.2) 

V = H%-\ (4.3) 

Now let and I7) be eigenstates of Hq with energies Ep and in the configuration 
of 3d 5 + 3d & L and in the 3d 6 configuration, respectively. These energies are defined from the 
ground state energy. The excited states containing a pair of a Ap electron and a ls-core hole, 
which is created by the El transition, may be given by p^aS^fi) . Also the excited states 
caused by the E2 transition may be given by |c) = s^It)- States |/3)'s span the space of 2352 
dimensions, and |7)'s span the space of 210 dimensions. Noting that the Is hole and the 4p 
electron are coupled to 3d — 2p electrons only through V, we have 

[G(z)]^ aY ee Miz-Hy^W) 
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- -l 

(z - E 1 + ei s )<5 7i y - V-y,vP G o(z - E^V^y , (4.4) 
where 

Vrfi* = (PlP^H^h), (4.5) 

G0{Z) = jk^z-e 4p M + e ls + i r> (46) 

with T being the core-hole life-time broadening width. The right hand side of eq. (4.4) means 
the inverse of the matrix whose components are written inside the brace. The inversion of the 
matrix is numerically carried out. 

Once we know [G(z)\ aiy(J1 i , we immediately obtain other components of the Green func- 
tion, 

[G(z)] avf 3, ai = (( 3 \4p' V aG(z)s (T \-f) 

= G (z-E p )J2v v( 3, Y [G(z)} aY , ai , (4.7) 

l' 

G(z) a7 , avf3 = ( 7 \siG(z)p^s a \(3) 

= Yt^WU^'V-f'^Goiz-Ep), (4.8) 

7' 

G(z) aT] p jaT] ,p, = {Plslp'^G^p'^SalP), 
= Go(z)5 v rf5[3 tl3 > 

+ Goiz-E^YVnpAGiz^ayVyrfpGoiz-Ep,). (4.9) 

77' 

The Green function is diagonal with the a variable. It should be noted here that the off- 
diagonal components, [G(z)] (TJ?j g )Cr7 and [G(z)] a ^ !art /3 given by eqs. (4.7) and (4.8), change their 
signs between the A\ and A-i sites, in accordance with the change of the effective coupling 
between Ap and 3d states. In eq. (4.9), the last term could not appear at the centrosymmetric 
sites, since it arises from the effective coupling between the Ap and 3d states which are not 
allowed within the same site. Note that, if a larger size of cluster is considered, the p-symmetric 
states could couple to such "3d" states through neighboring iron sites. 

Among many |/3)'s in the 3d 5 + 3d 6 L configuration, the lowest energy state \(3q) is taken 
into account in the following calculation. This may be justified when the presence of the pair of 
Ap electron and ls-core hole could not modify the 3d states through the Coulomb interaction. 
This observation simplifies greatly the analysis of the K-edge RXS in the next section. 
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5. X-ray absorption and scattering near the K-edge of Iron 

5.1 Transition matrix elements 

We need to consider two processes around the -ftT-edge; one is the El process that the 
Is core electron is excited to the Ap states, and the other is the E2 process that the ls- 
core electron is excited to the 3d states. These processes may be represented by transition 
operators, 

T™(j) = M lP ^s a , T™(j) = M 2 S ma s c , (5.1) 

where rj (= x, y, z) and m (= x 2 — y 2 , 3z 2 — r 2 , yz, zx, xy) are connected to the polarization of 
the incident photon. The annihilation and creation operators are defined with respect to the 
jth iron site. Since the Is state is well localized around the iron site, M\ and M 2 are evaluated 
by using atomic wavefunctions. We have 

Mi = iq J (^ p |a# ls )d 3 r 

i r°° 

= iq^ r 3 R 4p (r)R ls (r)dr = iAAQ x 1(T 3 , (5.2) 
v3 Jo 

M 2 = -q 2 j{^ z d x \{zx/2)\^ ls )d?r 
1 f'°° 

= -q 2 -j= \ r 4 R 3d (r)R ls (r)dr = -4.07 x 10" 4 , (5.3) 
zy 15 Jo 

where Ri 8 (r), i?3d(r), and i?4 P (r) are radial wavefunctions of the Is, 3d, and 4p states, respec- 
tively, which are calculated within the HF approximation. 32 We have inserted q ~ 3.6 x 10 8 
cm -1 for the x-ray wavenumber, which corresponds to the K-edge energy 7.12 keV. 

5.2 Absorption and XMCD spectra 

Although our main concern in this paper is the RXS spectra, we briefly discuss the ab- 
sorption spectra for looking over the whole -fT-edge region. 

The absorption coefficient may be given by a sum of contributions from each site, since the 
Is state is well localized at one atomic site. In general, it is decomposed into the contributions 
of the El-El, EI-E2, E2-EI, and E2-E2 processes: 

4» = Y,^K 1 Hm{n\T^\j)\g)5{u-E n + E g ), (5.4) 
= ^9\T^\j)\n){n\T^{j)W^-E n + E g ), (5.5) 

jn,(j 

= Y,^\T^Hj)\n)(n\T^\j)\9)^-E n + E g ), (5.6) 

A% m (u>) = Y.^\T^{j)\n){n\T^l (j)\g)S(u - E n + E g ), (5.7) 

where |g) and |n) represent the ground and excited states of the system with energy E g and 
E n . For example, when the x-ray is traveling along the z-direction with the polarization vector 
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along the x-direction, we need to set rj = x and m = zx. To include the life-time broadening 
of the core level, we replace the (5-function in eqs. (5.4)-(5.7) by the Lorentzian function with 
the full width of half maximum 2r = 1.6 eV. 

In the main K-e&ge region, the absorption coefficient is given by 

4£M = 2|Mx| 2 (-^j hnG (u), (5-8) 

where Gq(uj) is defined by eq. (4.6). It is expressed by the sum over k and can be replaced 
by the integral of the 4p DOS. It is known in many transition-metal compounds that the 
absorption spectra are well reproduced in the wide range 20 ~ 30 eV by means of the 4p DOS 
given by the band calculation. 35 In this paper, instead of carrying out the band calculation, 
we assume the 4p DOS rising from the energy corresponding to u = 7111 eV with the band 
width as large as 30 eV and sharp cutoff, so that it reproduces the experimental absorption 
spectra in the main K-edge region (see Fig. 3). 

Focusing on the contributions from the A sites in the pre-edge region, we have a more 
accurate form for A^L{u). Equation (5.8) is modified by including the last term of eq. (4.9), 

A^(u) = IMxI^Z^^H, (5.9) 

CT 

with 

^AWAjM = 2^ {\. G ( u )]*v'Po,<tvPo ~ [G'MW.WA.} • ( 5 - 10 ) 
Here G* is a complex conjugate of G. Only the lowest energy state |/3o) m the 3d 5 + 3d 6 L 
configuration, which is equivalent to \g), is considered by the reason explained in the previous 
section. The contribution of the last term of eq. (4.9) is, however, one order of magnitude 
smaller than that from the first term, and A^L{u) is practically determined by eq. (5.8). Note 
that AIUlo) could include the contribution of the 3d states at further neighbor iron sites if a 
larger cluster is considered. 

The contributions of the E1-E2 process, A^iu) and (u>), are canceled out after 
summing over the A sites, since they are proportional to [G^)]^^^ and [G(u)\ aijar) p at 
each site, and their signs change between at the Ai and A2 sites. Therefore, the breaking of 
centrosymmetry could not influence the absorption spectra. 26 The contribution of the E2-E2 
process is given by 

A^iui) = |M 2 | 2 Ysidld^Dll^i^ij'ldljg), (5.11) 
0-77' 

with 

D 22 W M = ^ {[G(u,)]* 7>7 - [G( W )]^}. (5.12) 

The upper panel in Fig. 3 shows the calculated absorption spectra in comparison with 
the experiment. 36 Any reliable theoretical estimates of the core-level energy are not available. 
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In addition, the K-edge energy is different for different experiments. 22-24 Therefore, we have 
tentatively adjusted the energy separation between the ls-core level and the Ap states. Since 
A^ m (oj) is limited within the pre-edge region, the spectra in the main i"C-edge region are 
dominated by A!A(lo). The band bottom of the 4p DOS corresponds to u = 7111 eV. A tail in 
A^(uj) due to r gives a substantial contribution in the pre-edge region, as shown in the inset 
in the figure. The total intensity in the pre-edge region is underestimated in comparison with 
the experiment. 36 



0.4 



Q 

o 



X 



0.000 



-0.001 

7100 




(b) XMCD 




7110 7120 
co[eV] 



Total 

E2-E2 

Exp. 

I I L_ 

7130 



Fig. 3. Absorption coefficient (upper) and XMCD spectra (lower) as a function of photon energy. The 
solid and dotted lines represent the calculated and the experimental spectra, 36 respectively. The 
inset is the decomposition of the total spectrum into A n (Lj) and A 22 (uj) in the pre-edge region. 



When the x-ray is traveling along the direction opposite to the magnetization, the absorp- 
tion coefficient is different between the right-hand and left-hand circular polarizations. The 
XMCD is defined by the difference between them. It is known that the XMCD is brought 
about by the SOI. We neglect the SOI on the 4p states, since its effect is expected to be very 
small in the pre-edge region. The lower panel in Fig. 3 shows the calculated XMCD spectra 
in comparison with the experiment. 36 The calculated difference is divided by the value at the 
peak of the main edge in the absorption coefficient. Since no scale is shown for the XMCD 
spectra in ref. 36, the experimental curve is drawn in arbitrary scale. The E2-E2 process gives 
the largest contribution. 

Note that these results are obtained for the A sites. For the B sites, the main i^-edge 
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spectra are the same, but the pre-edge spectra could be different. We need to consider the 
contribution from the B sites for quantitative comparison with the experiment. 
5.3 RXS spectra 

We consider the scattering geometry as illustrated in Fig. 4, where the incident x-ray with 
momentum k, energy u, polarization e is scattered into the state with momentum k', energy 
lo, polarization e'. We define the scattering vector by G = k' — k. 37 

scattered photon 




Fig. 4. Scattering geometry. Incident x-ray with momentum k and polarization a or 7r is scattered 
into the state with momentum k' and polarization a' or n'. 



By the same reason as the case of the absorption spectra, the RXS amplitude may be 
given by a sum of amplitudes from each iron site. Then the scattering amplitude per unit cell 
is expressed as 

F(G,u) = n>[i^(G)e- £ ' + X;i^41i( G . w )^ + E^' F "m( G ' w )^ 

+ E^'^VG.^+E^'^mCG.w)^], (5-13) 

m',r] m',m 

where the classical electron radius ro = e 2 / (mc 2 ) = 2.82 x 1CT 13 cm. The first term represents 
Thomson scattering, which may be estimated as 

F Th (G) = £ /°(G) «p(-iG • r,), (5.14) 

3 

where fj(G) is the atomic form factor with j running over not only iron sites but also oxygen 
sites. 

The remaining terms represent resonant scattering. They are defined by 

F c A /^(G,u,) =^// A Hc(exp(-iG • tj), A, A' = 1,2, (5.15) 

3 

where /^(lu)^^ is the resonant scattering amplitude at the jth iron site in the unit cell. For 
example, f} 2 {u) is defined by 

m,„. _ nu* t wmm^mi. (5 . 16) 
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Other components are similarly denned. In eq. (5.13), P% and Q m represent the geometrical 
factors of the El and E2 transitions for the incident x-ray with polarization = a or ir, while 
P% and Qm represent those for the scattered x-ray with polarization //' = a' or it'. Their 
general expressions are summarized in Appendix. 

The resonant terms at the A sites are expressed by using the resolvent given in Sec. Ill: 



/IV) 
/IV 



m'rj 



mC 2 \M 1 \ 2 Y J [G(^)]ar,',a V , 
a 

mc 2 M* l M 2 Y J [G(u;)U / ^( 1 \dl (7 \g), 
mc 2 M 2 *M 1 Y,{9\d m '*h)[G(Lo)U ,trr)l3 i 

(77 

mc 2 \M 2 \ 2 ^2(g\d m ' a \i)[G{uj)] (7 ^ cr y('y'\d ma \g). 



(5.17) 
(5.18) 
(5.19) 
(5.20) 



(777' 



Amplitudes fj^(oo) v ' m and fj^(u>) m ' v change their signs between the Ai and A 2 sites, respec- 
tively. They depend also on the direction of the local magnetic moment. A careful examination 
of eqs. (5.18) and (5.19) leads us to the expression, 
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(5.22) 



where (n x ,n y , n z ) represents the direction cosine of the local magnetic moment vector centered 
at each iron atom. Note that the local magnetic moment at the A sites is opposite to the 
total magnetization. The component a(u>), which is independent of the direction of the local 
magnetic moment, exists even in the absence of the SOI. On the other hand, b(oo), c(uj), and 
d(uj), which are one order of magnitude smaller than a(w), disappear without the SOI. All 
these components are appreciable only in a narrow pre-edge region. 
For the resonant terms at the B sites, f&~(u)) may be expressed by 

&\u)rfr, = mc 2 |M 1 | 2 2G M<W ( 5 - 23 ) 

where the contribution like the last term of eq. (4.9) does not exist. The f^(u>) and /^(lo) 
disappear because of centrosymmetry. We expect |/g 2 (u;)| << I/b 1 ^)! m the pre-edge region. 
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This contrast with the absorption coefficient, where the E2-E2 contribution is comparable 
to the El-El contribution. The reason is that the scattering amplitude is affected by Gq(lu) 
itself, whose real part is about two orders of magnitude larger than the imaginary part, while 
only the imaginary part contributes to the absorption coefficient. Note that, though it is 
small, /g 2 (w) could give rise to the magnetic scattering amplitude, which study is outside of 
the purpose of this paper. 

In the following, we analyze the RXS spectra at several Bragg spots, focusing on the a — a' 
channel. 

5.3.1 G = (002) and (006) 

For position vectors given by eq. (2.1), the phase factors exp(— iG ■ rj) are 1 at the Ai sites, 
— 1 at the A2 sites, and =fi,^fi,^fi,=fi,±i,±i,±i,±i,±i,±i,±i,±i,=fi,=fi,=fi,=fi,=fi (upper 
and lower signs correspond to (002) and (006), respectively) at the B sites. Thomson scat- 
tering amplitude as well as all the resonant terms are canceled out except for / A 2 (u;) and 
f^-(uj), due to the phase factors. Therefore, these Bragg spots are suitable to investigate non- 
centrosymmetric effects on the RXS. Several experiments of RXS have actually been carried 
out on these spots, 21 ' 22 but the dependence of the spectra on the magnetization direction has 
not been studied yet. We calculate the RXS spectra and analyze such dependence. 

We consider the situation that the scattering plane contains a vector (1,-1,0) (see Fig. 4). 
The geometrical factors P CT , P u \ Q a , and Q u ' are given by putting ip = 7r/4 in the expressions 
in Appendix. We assume that the local magnetic moment on the A sites is along the n = 
(n x ,n y ,0) direction. Then, using eqs. (5.21) and (5.22), we have the scattering amplitude in 
the a — a 1 channel, 

F(G,U) = r ^ [^'^m( G > u )Qm + Qm F mr)( G i U ) P % 
r},m 

= 8r [2 sin 6a(uj) + V2(n x - n y ) cos 6b(u)}, (5.24) 

where Bragg angle 9 is 12.0 and 38.5 degrees for (002) and (006), respectively. Let I + (G,u) 
and I-(G,uj) be the intensities per unit cell for the direction of the magnetic moment n and 
the reverse, respectively. Then the average and the difference of the intensities are given by 

r- CT '(G,u,) = i(/r CT '(G,u,) +J p- CT '(G,u,)) 

= 256r 2 {sin 2 e\a(uj)\ 2 + ^cos 2 6(n x -n y ) 2 \b(u;)\ 2 }, (5.25) 
Ar- a '(G,uj) ee I°- a '(G,uj) -/r CT '(G,u;) 

= 128v / 2r 2 sin26»(n :r - n y )[a(uj)*b(uj) + c.c.]. (5.26) 

Since b(uo) is one order of magnitude smaller than a(cj), the average intensity I <T ~ a ' (u) is 
dominated by |a(a;)| 2 . On the other hand, the difference spectra arise from the interference 
between the terms of a{uj) and b(uj). 
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Figure 5 shows I a ~ a ' (G,uj) and AI a - a ' (G, u) calculated with n = (1/^2,-1/^/2,0). 
This magnetization direction corresponds to the magnetic field applied along the [—1,1,0] 
direction. The I a ~ a '(G,uj) is concentrated in a narrow pre-edge region, and becomes larger 
for G = (006), in consistent with the experiments. 21 ' 22 The AI a ~ a ' (G, to) is relatively large, 
only one order of magnitude smaller than the average intensity. This intensity difference is a 
consequence of breaking both the local centrosymmetry and the time-reversal symmetry. It 
would not be hard to detect such a difference. 
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Fig. 5. Scattering intensity as a function of photon energy in the pre-edge region in the a— a' channel. 
Panel (a) shows the average intensity. Panel (b) shows the intensity difference when the magnetic 
field is applied along the [—1, 1, 0] direction and the reverse, which is divided by the peak value 
of the average intensity. The solid and broken lines represent the intensities for G = (002) and 
(006), respectively. The inset on panel (a) represents the experimental curves taken from ref. 22. 



5.3.2 G = (222) 

For position vectors given by eq. (2.1), the phase factors exp(— iG ■ rj) are 1 at the Ai 
sites, —1 at the A2 sites, and i at all B sites. Therefore this spot is not prohibited. 

Thomson scattering gives the largest contribution; the contribution from the A sites are 
canceled out, but those from the B sites and oxygen sites remain, resulting in 

F Th (G) = lGifa - 32i/o, (5.27) 
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Table II. F Th (G) and F 11 (G,oj ) with oj being the pre-edge absorption peak. 



hkl F Th (G) [F n (G,co )} m 
(222) i 123.1 0.68 -i 10.24 

(333) 148.7 -i 148.7 -32.52 + i 5.78 



where /b = (l/2)[/(Fe 2+ ) + /(Fe 3+ )] is the form factor of iron at the B sites and fo = / (0 2 ~) 
is that of oxygen. They are evaluated from the atomic values tabulated in ref. 38. As regards 
the resonant terms, fj^(cu) and fj^(uj) are canceled out at the A sites due to the phase factor, 
but the contribution from the B sites remains. We have 

[f U (G, W )]^ = F U (G, W )^ = 16i/ B V)<W ( 5 - 28 ) 

where f^-(u}) is evaluated from eq. (5.23). Note that the atomic calculation of /^(w) would 
contain large errors in the pre-edge region, because a single level of 4p states in an atom 
changes into an energy band with width as large as ~ 20 eV in solids. We list the calculated 
values in Table II. The Thomson scattering amplitude is much larger than the resonant term. 
These values are much larger than those reported in ref. 24. Thus, the average intensity is 
given by 

/ CT - CT '(G,^) =r 2 |F Th (G)+F n (G,^| 2 . (5.29) 

Next we analyze the dependence on the direction of applied magnetic field in accordance 
with the experiment. 24 The scattering plane is set to contain a vector (1, —1, 0) with applying 
magnetic field along the [112] direction and the reverse. The geometrical factors P u , P u ' , 
Q a , and Q a are given by putting ip = 7r/4 in the expressions in Appendix. Substituting 
(Tl/V6,Tl/v / 6,±2/V6) for (n x ,n y ,n z ) in eqs. (5.21) and (5.22), we obtain the scattering 
amplitude in the a — a 1 channel, 

[< [^ 12 (G, u)] vm Q° n + Qt [F 21 (G, u)] mv P%] = 8 cos e[5 a(uj) ± 6 lC (u) ± S 2 d(u)}, (5.30) 

i),m 

where the upper (lower) signs correspond to upper (lower) signs of (n x ,n y ,n z ), and 5q = 
-(2/3)(l/ v / 2 + 2/V3), Si = -V2/3, 5 2 = (2/3)(l - ^/2/3). Bragg angle 9 is 21.1 degrees. 
Since the direction of the local magnetic moment on the A sites is opposite to the direction 
of the applied magnetic field, we define the intensity difference as the value with the upper 
sign for n minus the value with the lower sign for n. As a consequence, we have 

Af" ff '(w) = 16rgcos0{[F Th (G)* + F^ 1 (lo)*][5ic{uj) + 5 2 d(uj)} + c.c.}. (5.31) 

The intensity difference arises from the interference between the term of F Th (G) + Fq 1 (G, uo) 
and the terms of c(u>) and d(u). 

Figure 6 shows the relative intensity difference AI a ~ a> (G,Lu)/I a ~ a ' (G,lj) in the pre- 
edge region. The spectral shape takes a peculiar "dispersion" form. This may be explained 
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as follows. The factor F Th (G) + F 11 (G,u>) at (222) spot is very close to a pure imaginary 
number as shown in Table II. Another factor c{uj) or d(uj) is given by a resolvent matrix 
element multiplied by M\M2 which is a pure imaginary number. Thereby the product of the 
two factors in eq. (5.31) becomes a resolvent matrix element multiplied by a real number. By 
adding its complex conjugate, eq. (5.31) becomes proportional to the real part of the resolvent 
matrix element. The real part of the Green function usually takes a "dispersion" form as a 
function of energy. 

The intensity difference in the experiment by Matsubara et al 2A is different from the one 
calculated here. It extends over a region much wider than the region of the pre-edge absorp- 
tion spectra with an order of magnitude larger intensity, which behavior is quite unusual. 
Also the shape is different from a "dispersion" form. Matsubara et al. claimed that the spec- 
tra they found arise from a "magnetoelectric" amplitude, that is, a consequence of breaking 
both the centrosymmetry and the time-reversal symmetry. According to the present analysis, 
this claim has no ground. In this connection, we would like to draw attention to the simi- 
lar RXS experiment for MnCr204, where the shape and strength quite similar to the curve 
calculated above have been observed at the Mn pre-K-edge. 27 In this material, Mn atoms 
occupy at the A sites in spinel structure. Since the Mn pre-edge spectrum selects only the A 
site contribution, this experimental result indicates that the calculated spectra correspond to 
the "magnetoelectric" signal. Note that a similar "dispersion" form of the spectra has been 
observed 16 and theoretically analyzed 17 ' 18 at the Fe K edge of GaFeO"3. 

Finally we comment on what happens on the intensity difference when the scattering 
vector is reversed. Different from the conventional case, the signal is reversed with keeping 
the shape, as shown in Fig. 6. The F l2 (G,uj) and F 12 (G,lo) are unaltered because the phase 
factors at the A sites are the same with reversing G. On the other hand, -F Th (G) + -F 11 (G, lo) 
changes its sign, because the phase factors at the B sites is changed from i to —i, resulting in 
the sign change in eq. (5.31). 
5.4 G = (333) 

The phase factors exp(— iG ■ rj) are 1 at the Ai sites, — i at the A2 sites, 
exp(i37r/4), exp(i37r/4), exp(i37r/4), exp(i37r/4), exp(— in/ 4), exp(— in/4), exp(— in/4), 
exp(— in/4), exp(— in/4), exp(—in/4), exp(—in/4), exp(—in/4), exp(in/4), exp(in/4), 
exp(in/4), exp(in/4) at the B sites, respectively. The Thomson scattering amplitude and 
the resonant term i ?11 (G,oj) are given by 

F Th (G) = 4(l-i)(/ A + V2fe), (5.32) 

F n (G,a;V„ = F n (G, u) <V„ = 4(1 - i)[f{\u;) + >/2fg(u)] (5.33) 

These amplitudes are evaluated by using atomic form factors for /a and /b and eq. (5.23) for 
fj^(io) and /^(lo). The results which are listed in Table II. In addition, we have the E1-E2 
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Fig. 6. Relative intensity difference AI a a (G,Lo)/I a a (G,u>) as a function of photon energy at 
G = (222) and G = (222) in the pre-edge region. 



term, 

£ [rf [F 12 (G, Lo)} vm Q^ + Qi\F 21 (G, u,)\ mv P°\ = 4(1 + i) cos 0[6 a(uj) ± 8 lC {u) ± <5 2 d(o;)], 

(5.34) 

with Bragg angle 6* = 32.6 degrees. As a result, we obtain the intensity difference as 

AI a - a '(io) = 8r% cos0{[F Th (G)* + F^(uj)*}(1 + i)[$ic(w) + <5 2 d(w)] + c.c.}. (5.35) 

Since F Tll (G)* is proportional to the right hand side of eq. (5.35) is nearly proportional 

to the real part of a resolvent matrix element. Therefore we would expect a "dispersion" form 
of spectra as a function of photon energy. 



0.0001 



3, 

1- 0.0000 

o 

< 



-0.0001 



7100 



t 

1 \ 
1 \ 
1 1 
/ 1 

; 1 

/ 

/ 


K. 


: 1 


r 

I 

1 / 
1 / 
1 f 
1 ( 

1/ G=(3,3,3) 

v -- - G=(-3-3-3) ^ 



7110 
co[eV] 



7120 



Fig. 7. Relative intensity difference AI a a (G,Lo)/I a a (G,w) as a function of photon energy at 
G = (333) and G = (333) in the pre-edge region. 



Figure 7 shows the relative intensity difference AI a CT '(G, uj)/I cf a ' (G,u) thus evaluated. 



20/26 



J. Phys. Soc. Jpn. Full Paper 

The calculated value has the same size of magnitude as the experimental one 24 and is nearly 
half of the calculated one for (222). Note that the experimental value at (222) is about two 
order of magnitude larger than the values at (333). 
5.5 G = (444) 

The phase factors exp(— iG ■ rj) are 1 at the A sites and —1 at the B sites. Therefore, the 
E1-E2 terms -F 12 (G, uj) and F 21 (G, u>) vanish due to the cancellation between the Ai and A2 
sites. In the experiment, 24 the intensity dependence is found negligible. If we take seriously 
this fact, it means that the contribution of the magnetic scattering amplitude F 22 (G,uj) is 
quite small. 

6. Concluding Remarks 

We have studied how the breaking of centrosymmetry affects the RXS spectra through a 
microscopic calculation for magnetite. The centrosymmetry is locally broken at tetrahedral 
(A) sites. In such a circumstance, the 4p states strongly hybridize with the 3d states through 
neighboring oxygen 2p states, giving rise to the non-vanishing contribution of the E1-E2 
process in the RXS spectra. This observation is substantiated by introducing a microscopic 
model of a FeC>4 cluster with the Ap states forming a band and the 3d states forming multiplet 
structures. We have calculated the RXS spectra with the help of the resolvent formalism. It 
is shown that the hybridization changes its sign between the Ai and A2 sites and accordingly 
the local amplitude from the E1-E2 process changes its sign. This sign change causes non- 
vanishing RXS intensities at the forbidden spots (002) and (006). The spectra are concentrated 
in a narrow pre-edge region with intensities larger at (006) than at (002), in agreement with 
the experiment. In addition, we have carefully analyzed the scattering matrix for the E1-E2 
process, which depends on the direction of the applied magnetic field. Such dependence is 
only possible when both centrosymmetry and time-reversal symmetry are broken. Through 
this analysis, we have obtained large dependences of intensity at (002) and (006) spots. We 
hope that this dependence could be observed in future experiments. 

We have also analyzed the dependence on the direction of the applied magnetic field at 
(222), (333) and (444) spots in connection with the experiment. These spots are allowed 
with large Thomson scattering amplitudes. Having calculated the intensities for two opposite 
directions of the applied magnetic field, we have obtained their difference with the same 
order of magnitude at both (222) and (333) spots but negligible difference at (444). The 
intensity difference is found to has a "dispersion" form as a function of photon energy, which 
is concentrated in a narrow pre-edge region. In the experiment by Matsubara et a/., 24 however, 
the intensity difference at the (222) spot is distributed over the region much wider than the pre- 
edge region with the spectral shape quite different from the "dispersion" form. The observed 
intensity difference at (222) is two orders of magnitude larger than the one at (333). These 
behaviors seem unusual and hard to explain. Matsubara et al. claimed that the difference at 
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the (222) spot arises from a "magnetoelectric" amplitude, that is, a consequence of breaking 
both centrosymmetry and time-reversal symmetry. This claim has no ground, according to 
the analysis in this paper. A "dispersion" form of spectral shape has been observed in the 
experiment at the Mn pre-i^T-edge in MnCr204. 27 Since Mn atoms are occupying at the A 
sites in spinel structure, the pre-if-edge signal selects the contribution from the A sites. This 
experiment suggests that the calculated spectra correspond to the "magnetoelectric" signal 
in the magnetite. The large signal at the (222) spot might be related to B sites. We hope 
experiments in future clarify the situation in magnetite. 
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Appendix: Geometrical factors 

We briefly describe a derivation of geometrical factors. A geometrical setting of x-ray 
scattering adopted in the present work is shown in Fig. 4. We introduce three coordinate 
systems, (x',y',z'), (x" ,y" , z"), and (x'" ,y"' , z'"). The first coordinate system is attached to 
the incident (scattered) photon with its z' axis being parallel to k(k'). Its x' and y' axes are 
perpendicular and parallel to the scattering plane, respectively. The second coordinate system 
is used as the definition of the origin of the azimuthal angle tp. Its z" axis is aligned to G 
direction and y" axis is in the scattering plane at ip = with x" = x'. The third coordinate 
system is fixed to the crystal. 

These three coordinate systems can be connected by the Euler rotation with the choices 
of appropriate Euler angles. Here we use the same definition of the Euler rotation adopted in 
Rose's book. 39 From (x" , y", z") to (x',y',z') coordinate systems, the Euler angles are chosen 
as ^ ± 9, — ^) where 9 represents the Bragg angle. The upper (lower) sign is for the incident 
(scattered) photon. Hereafter, we restrict our discussion on the incident photon case alone, 
since the results for the scattered photon are obtained by replacing every 9 with —9. The Euler 
angles of the rotation from the (x" ,y" , z") coordinate system to the (x'", y'" , z 1 ") coordinate 
system are given by (a, (3, 0) where a and (3 are the azimuthal and the polar angles of G, 
respectively. 

In order to calculate the geometrical factors, we start with writing down the basis corre- 
sponding to the a and it polarizations. For dipole transition, it is simple since e a = e x r and 
e 71 " = —e y > hold in the present setting. Here ej denotes the unit vector directed to j axis. 
Then, the geometrical factors {Pj 1 } = 7r) are defined by the following relations. 



e - . r = x' = P a x x" + P°y"' + P a z z"', 
e n ■ r = -y' = P^x'" + P?j/" + P?z' : 



(A-l) 



(A-2) 
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where r is an arbitrary position vector. Similarly, the geometrical factors {Qn} are defined by 
the quantity (k ■ r)(e ■ r) appearing in the multipole expansion of the scattering amplitude. 
Here n = 1, 2, 3, 4, and 5 correspond to the quadrupole basis x 2 — y 2 , 3z 2 — r 2 ,yz, zx, and xy, 
respectively. By noticing the fact that k = |k|e 2 /, we define {Qn} in the following relations. 



(kT)(6*T) oc 4 = 

71=1 

5 

(k • r)(e 7r • r) oc -z' 3 = ^ Qn^'n , 



(A-3) 
(A-4) 



n=l 



where z\ = ^(x 2 — y 2 ) Z2 = \{3z 2 — r 2 ), z 3 = V3yz, z^ = \/3zx, and z§ = V3xy. 

From eqs. (A-l) ~ (A-4), the geometrical factors are evaluated by expressing x', y' in terms 
of x'", y'", z'" for the dipole transition and z ; 4 , z' 3 in terms of z" n ' for the quadrupole transition, 
respectively, with the help of rotation matrix. The final results for the a polarization with the 
incident photon are as follows: 

P£ = cos a cos (3 cos ip + sin a sin ip, 

sin a cos (3 cos ip — cos a sin ip, 



y 



P° = — sin/3cos'0, 



1 



(A-5) 
(A-6) 
(A-7) 



Q1 = - sin 6* (— cos V>cos 2a sin 2/3 — 2 sin ?/> sin 2a sin /?) 



Ql 



+ - cos [sin 2ip cos 2a (1 + cos 2 (3) — 2 cos 2ip sin 2a cos (3] 
= sin # cos ^ sin 2/3 H — — cos 6 sin 2^ sin 2 (3, 



Q 3 = sin#(— cos ■;/; sin a cos 2/? + sin ?/> cos a cos /3) 

+ - cosO(— sin 2ip sin a sin 2(3 — 2 cos 2-0 cos a sin /3) , 

Q4 = sin#(— cos V 7 cos a cos 2/3 — sin f/' sin a cos /3) 

+ i cos6 l (— sin 2^ cos a sin 2/3 + 2 cos 2-0 sin a sin /3) , 

Q5 = - sin 6 (— cos ip sin 2a sin 2(3 + 2 sin ?/> cos 2a sin (3) 

+ i cos 6 1 [sin 2^ sin 2a (1 + cos 2 /3) + 2 cos 2-0 cos 2a cos 0\ 



(A-8) 
(A-9) 

(A-10) 

(A-ll) 

(A-12) 



For (004£ + 2), 


putting c 


x = (3 = 0, 


we 


sin 2^ cos 9, Q^ 


= 0,Qf = 


- sin ip sin # 


, Q4 


a = 7r/4 and (3 


= sin -1 y 


2/3, we have 




- 1 X 


- /I 


cos 




pa 

y 


- /I 


cos 



(A-13) 
(A-14) 
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Qi 



•2 

'1 
1 fl 



cos V', 

— cos 2ip cos 9 — \[2 sin ip sin 9 , 
sin 2ip cos 9 + \pl cos -0 sin 9 , 
[— \/6 cos 2ip cos — y/2 sin 2-0 cos 



3 

+ cos ip sin sin V' sin 0] , 

1 fl 

3 

+ cos ^ sin sin ^ sin 0] , 



= - J- [y/E cos 2ip cos - >/2siii2^cos0 

o V ^ 



Q5 = -— - [\/2 sin 2?/> cos 6* _ cos sin ( 



Full Paper 
(A-15) 
(A-16) 
(A-17) 

(A-18) 

(A-19) 
(A-20) 



The expressions for the scattered photon are obtained by replacing 9 with —9 in the above 
expressions. 
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